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CONPAfilSON OF ANALYTICAL AND EXPERIMENTAL SUBSONIC STEADY- AND UNSTEADY- 
PRESSURE DISTRIBUTIONS FOR A HIGH-ASPECT-RATIO SUPERCRITICAL WING 
MODEL WITH OSCILLATING CONTROL SURFACES 
By HI Ilian E* McCain 

SUMMARY 

The results of a comparative stuc(/ using the unsteac(y aeroc^ynanlc lifting- 
surface theory, knoMn as the Doublet Lattice method, and experimental subsonic 
steady- and unsteady-pressure measurements, are presented for a high-aspect- 
ratio supercritical wing model. Comparisons of pressure distributions due to 
w;.'.g angle of attack and control -surface deflections were made. In general, 
good correlation existed between experimental and theoretical data over i>x)St of 
the wing planform. The more significant deviations found between experimental 
and theoretical data were in the vicinity of control surfaces for both static 
and oscillatory control -surface deflections. 

INTRODUCTION 

The application of active controls technology, especially in conjunction 
with such advanced aerodynamic features as winglets and supercritical airfoils, 
can provide improved aircraft characteristics with relaxed static stability, 
load alleviation, and flutter suppression systems. The design and analysis of 
active control systems for energy-efficient transports frequently require multi- 
purpose computer programs that include unsteady aerodynamic calculations. 
Programs currently used at Langley Research Center for the synthesis and 
analysis of active controls include !SAC (ref. 1), SYNPAC (ref. 2), DYLOFLEX 


(refs. 3 and 4), and the aerodynamic energy method (ref. S), all of Mhich use 
the Doublet Lattice method for unsteady aerodynamics. An active con'.rol system 
design Is very sensitive to the quality of the aerodynamic results eii|>1cyed. 

To investigate the effect" of oscillating control surfaces on unsteady 
aerodynamics, a series of wind-tunnel tests were conducted by the National 
Aeronautics and Space Administration (NASA) In the Langl^ Transonic Dynamics 
Tunnel (TDT), using a semispan model of a high-aspect-ratio supercritical wing. 
The primary objective of these tests was to provide a comprehensive data base of 
measured transonic unsteady pressures for a wing representative of those 
envisioned for an energy-eff icient transport. A secondary objective was to 
provide similar data at subsonic conditions. The results from these tests are 
contained In references 6 and 7. These results are In the form of both steady 
and unsteady pressure distributions for variations In Mach number, Reynolds 
number, wing angle of <ttack, control -surface deflection angle (static and 
oscillating), and control -surface oscillation amplitude and frequency. 

The purpose of this paper is to compare the experimental pressure 
distributions available at Mach number 0.60 with pressure distributions 
calculated by the Doublet Lattice method. The Doublet Lattice method was chosen 
because of its extensive use In active controls synthesis and analysis. 

Lxclusive use of Deublet Lattice results In this paper should not be Interpreted 
as disregarding other available methods, such as the Kernel function method of 
RHOIV (refs. 8 and 9), or the more advanced methods of small -perturbation 
or full-potential equation fonmiatlons (reviewed in refs. 10, 11, and 12). 

Since it is well known that Doublet Lattice, as well as other Hftlng- 
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surface theories* overpredicts control -surface aerodynaulcs* correction factors 
are frequently applied to reduce this error (ref. 13). It Is hoped that the 
results of this paper Mill assist In the development of empirical correction 
methods that can be applied to the Doublet Lattice calculations. 

SYMBOLS 


AR aspect ratio, bo^/S 

bo/2 wing semi span, m 

b wing root semi chord, m 

c streamwise local chord, m 

Cav wing average chord, m 

c^ section lift coefficient 

c^i section pitching-moment coefficient about the leading edge 

Cp pressure coefficient 

C'p lifting-surface steady-pressure coefficient, Cp^ ^ - 

^Pu.s. 

AC|^/Aa Increment In section lift coefficient per change In angle of 

attack, deg"^ 

aC|^/a 6 Increment In section lift coefficient per change In control- 

surface deflection, deg"^ 

ACm/Aa Increment in section pitching-moment coefficient about the 

leading edge per change In angle of attack, deg"l 

AC'p/Aa increment In lifting-surface pressure coefficient per change 

in angle of attack, deg"l 

AC'p/A6 Increment In lifting- surf ace pressure coefficient per change 

In control -surface deflection, deg“^ 
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magnitude of lifting-surface unsteac(y-pre$sure coefficient 

frequency of oscillating control surface, Hz 

reduced frequency, bu/V 

free-stream Mach number 

free-stream dynamic pressure, kPa 

Reynolds number based on wing average chord 

total wing planform area, 

thickness-to-chord ratio 

free-stream velocity, m/sec 

streamwise coordinate, m 

fraction of local streamwise chord 

spanwise coordinate, m 

vertical coordinate, positive up, m 

wing angle of attack, deg 

control -surface deflection angle, positive trailing edge 
down, deg 

change in wing angle of attack, deg 

change in control -surface deflection angle, deg 

leading-edge sweepback angle, deg 

fraction of semi span, Zy/bg 

phase angle of unsteady pressure, referenced to control - 
surface motion (negative for pressure changes lagging 
the control -surface motion) 

circular frequency of oscillating control surface, rad/sec 
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Subscripts and Abbreviations 


c.s. 

control surface 

l.e. 

leading edge 

l.s. 

lower surface 

ref 

reference 

u.s. 

upper surface 


MODELS 

Wind-Tunnel Model 

A sketch of the wind-tunnel model Is presented In figure 1. The model has 
an aspect ratio of 10.76, a leading edge sweepback angle of 28.8®, and a 
semi span of 2.286 meters. total of 252 static- pressure orifices were 
installed on the wing, with half on the upper surface and the other half on the 
lower surface. The table in figure 1 lists the designated spanwise stations for 
each of the nine chordwise sets of orifices. A total of 164 dynamic-pressure 
transducers were installed on the wing at closely-corresponding locations to the 
static-pressure orifices. The model was equipped with a total of 10 oscillating 
control surfaces, 5 at the leading edge with hinge lines on the 15 percent chord 
and 5 at the trailing edge with hinge lines on the 80 percent chord. Only 2 of 
the 10 control surfaces were considered in this study: an inboard trailing edge 

control surface and an outboard trailing edge control surface (identified as 
nos. 6 and 9 in refs. 6 and 7). The inboard control surface is located between 
10 and 24 percent semi span, and the outboard control surface is loacted between 
59 and 79 percent semi span. The cross-sectional shape of the model consists of 
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NASA supercritical -airfoil sections of varying chord length and thickness as 
shown In figure 2. Further details of the wind-tunnel model. Including airfoil 
shape quality and planform rigidity are described In references 6, 7, and 14. 

Analytical Model 

An aerodynamic model was generated for use In the subsonic unsteady 
lifting-surface theory known as the Doublet Lattice method (ref. 15). The 
arrangement of aerodynamic boxes representing the wind-tunnel model Is shown In 
figure 3. To provide more calculated pressure points for comparative purposes, 
the chordwise and spanwise distribution of aerodynamic boxes were Increased over 
the planform areas near and on the control surfaces. The aerodynamic boxes In 
figure 3 with asterisks identify the control surfaces and the cross-hatched 
strips Identify the locations corresponding to the nine semi span stations shown 
In figure 1. There were 42 streamwise strips, a total of 325 aerodyO'^mlc boxes, 
used to comprise the model layout. Ii creating this box layout, an attempt was 
made to keep the aspect ratio of each box as close to 1.0 as possible. 

EXPERIMENTAL DATA 

All the experimental data presented herein were obtained at the test 
conditions of M « 0.60, R ■ 2.2 x 10^ (based on the wing average chord), and 
q - 3.0 kPa. The tests Included two angles of attack of zero degrees and 
approximately 2.85 degrees (corresponding to the cruise condition). 

Steady-Pressure Data 

A sumnary of the conditions at which the steady-pressure data were taken Is 
presented in Table I. The test point numbers used to identify the data In 


6 



reference 7 are listed In the table for convenience. Test conditions for both 
Inboard and outboard control surfaces Included seven control-surface deflection 
angles ranging from -6* to 6*. 

The results of reference 14 Indicated that the wind-tunnel model was 
essentially rigid and that pressure-measurement results were not significantly 
Influenced by model flexibility. Therefore, only four rigid body modes (plunge, 
pitch. Inboard, and outboard control -surface deflections) of the model were 
Included In the analysis. Reduced frequencies of 0, 0.14, 0.27, and 0.41 were 
chosen, corresponding to those at which experimental data were available In 
reference 7. 

The experimental data from ref. 7 consisted of tabulated steady-pressure 
coefficients, Cp, on both upper and lower surfaces of the wing with the 
corresponding lifting-surface steady-pressure coefficients, C'p « ^pi.s. " 

Cpu 5 . The section lift coefficient, c^, and section pitching-moment 
coefficient about the leading edge, c,n> spanwise stations 

were obtained by numerically Integrating equations (1) and (2), respectively. 

0 ) 

Cm C'p X dx (2:- 

c 


Calculations were also made for the Incremental changes In these coefficients 
due to angle of attack or control -surface deflection changes, as follows: 
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(5) 


( 6 ) 


(7) 


Typically, the referenced quantities corresponded to the zero-valued test 
conditions. 


Unsteady-Pressure Data 

A summary of the conditions at which the unsteady-pressure data were taken 
is presented in Table II. Again, the test point numbers identify the data 
obtained from ref. 7. Each control surface was oscillated about a zero mean 
deflection angle with an amplitude of ±4® for three frequencies of oscillation 
(5, 10, 15 Hz). Depending on the exact tunnel speed for a given test point, 
the corresponding reduced frequencies varied slightly about the average values 
of 0.14, 0.27, and 0.41, respectively. The unsteacly-pressure measurements are 
presented in the form of the magnitude of the lifting-surface unsteaciy-pressure 
coefficient, jC'pj. and phase angle, All phase angles were referenced to 
the control-surface motion, with negative values for pressure changes lagging 
the control -surface motion. Although unsteacly-pressure measurements were made 
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at all nine semlspan stations, only the chordwise distributions at two locations 
(one near the center of each control surface— rows 1 and 6 In figure 1) were 
considered In this study. 


ANALYSIS 

The Doublet Lattice formulation solves the linearized acceleration or 
pressure potential -flow equations on zero thickness lifting surfaces at subsonic 
speeds with nonplanar boundary conditions. The Doublet Lattice method (ref. 15) 
was used to generate the theoretical steady end unsteady aerodynamics herein. 

The calculations were performed by the version of the Doublet Lattice program 
which Is used in a NASA computer program system known as ISAC (Interaction of 
Structures, Aerodynamics, and Controls, ref. 1). As stated previously, four 
rigid-body modes (plunge, pitch. Inboard, and outboard control -surface 
deflections) were used at four reduced frequencies at a Mach number of 0.60. 

For each mode and at each reduced frequency, the output from the Doublet 
Lattice program consists of complex lifting-surface pressure coefficients on 
each aerodynamic box. Since the program performs the necessary numerical 
Integrations internally, the complex section lift and moment coefficients are 
also listed. At zero reduced frequency (steady), the Imaginary parts of these 
complex quantities are zero. The real and imaginary parts of tht unsteady 
quantities were converted to magnuudes and phase angles for direct comparison 
to the experimental values from ref. 7. 

COMPARISON OF ANALYSIS AND EXPERIMENT 
Steady-Pressure Results 

Comparisons of the incremental lifting-surface pressure distributions for 
an Incremental angle of attack are shown in figures 4(a) to 4(1) at each of the 
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nine setnispdn st.dt.ions. The two sets of ex|>eriiiient<ii ddtn correspond to the 
zero-deflect lon-tinyle data taken for both inboard and oiitboard control surfdces 
at an incremental arujie of attack. At each semispan station and for 
approximately the forward 4U percent of the local chord, tlie Doublet Lattice 
data noticeably underpredict the maynitude of Che chordwise pressure 
distribution. addition, for about the aft 40 percent of the local chord, it 
slightly overpredicts the magnitude. The deviations shown between the 
experimental and theoretical data are typical of airfoil thickness effects 
(viscosity), ref. 8. 

The next comparison illustrates the integrated result of the local 
incremental lifting-surface pressure distributions for section lift and 
pitching-moment coefficients, shown in figures 5(a) and 5(b). Figure 6 presents 
the spanwise distribution of local aerodynamic center locations. In figures 5 
and 6, the two sets of experimental data show good agreement between each other 
for measurements from approximately 25 to 80 percent semi span. However, there 
is limited scatter in the experit.iental data inboard of 25 percent semispan and 
outboard of 80 percent semispan, possibly due to fuselage-body and tip effects, 
respecti vely. 

The incremental lifting-surface pressure distributions for incremental 
cor'crol -surface deflections are shown in figures 7 and 8. Figure 7 presents 
comparisons for the inboard control -surface data at semi span station n ■ 0.19; 
figure 8 presents comparisons for the outboard control surface data at semi span 
station >i => 0.71. Each figure contains data at the two angles of attack a’'d 
for three positive and negative incremental control -surface deflection angles 
(Zi6 “ ±2", t4®, ±6®). The magnitudes of the experimental pressure coefficient 
due to positive (trailing edge down) control -surface deflections are less than 
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those due to equivalent negative deflections. Although experimental chordwise 
pressure distributions for either positive or negative control -surface 
deflections show similar trends, there Is better correlation with the theory for 
the experimental negative deflection data. For both control surfaces, the 
largest deviation between experimental and theoretical data occurs aft of the 
hinge line. The more significant difference occurs on the outboard control 
surface. Note that the experimental data for positive outboard control -surface 
deflections exhibit a rise In pressure at the 95 percent chord location after 
dropping to a minimum pressure at the 90 percent chord location. This effect Is 
possibly due to Increased spanwise flow or separated flow toward the trailing 
edge of the deflected control surface. The comparison between experimental and 
theoretical data forward of each control -surface hinge line shows good 
correlation in trend, but not In magnitude. The Doublet Lattice data in this 
region of the chordwise pressure distribution usually underpredicted the inboard 
measurements and slightly overpredicted most of the outboard measurements. This 
particular deviation being a percentage difference in magnitudes could be 
adjusted by correction factors to the Doublet Lattice aerodynamics, similar to 
those used on control surfaces, ref. 13. There were relatively snail 
differences overall between the respective data sets at each angle of attack. 

The increjnental section lift coefficient distributions for the incremental 
control -surface deflections are shown in figure 9. The averaged experimental 
values shown are for the integrated results of the local incremental lifting- 
surface pressure distributions due to the six incremental control -surface 
deflections (a 6 ■ tZ®, t4®, »6“). The experimental scatter is indicated by the 
vertical lines through each symbol. All the experimental incremental section 
lift coefficients shown are based on the zero angle of attack data. The 
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experimental averages for the Inboard control -surface data shM a good 
con^)ar1son with the theoretical results. The experimentat averages for the 
outboard control surface do not compare well with theory. As stated previously, 
the possible effect of spanwise flow or separated flow in this region could be 
significant. In addition, from the data presented herein and other data 
contained in references 6 and 7, there is evidence that small angle of attack or 
control -surface deflection changes can produce large pressure changes at the 
outboard semispan stations, n * 0.71, 0.78, 0.81, and 0.92. 

Unsteady-Pressure Results 

Comparisons of the chordwise unsteady-pressure distributions, in the form 
of magnitudes and phase angles, are presented in figures 10 and 11. For the 
three oscillating frequencies of 5, 10, and 15 Hz, the magnitudes and phase 
angles of the lifting- surface unsteady-pressure coefficients are compared at 
semispan station n = 0.18 (figures 10(a)-(c)) for the inboard contr <1 -surface 
data, and at semispan station n * 0.71 (figures 11(a)- (c)) for the outboard 
control -surface data. The experimental data for both angles of attack (zero and 
cruise) exhibited only minor scatter. Similar to the steady-p<"essure 
conparisons, the Doublet Lattice data overpredicted the experimental unsteady- 
pressure magnitudes aft of the control -surface hinge lines, with the more 
significant deviation on the outboard control surface. In addition, there is an 
overprediction of the experimental data at the 75 percent chord of the outboard 
control surface (immediately forward of the hinge line). The two trailing edge 
(x/c * 0.90, 0.95) experimental unsteady -pressure measurements also exhibited 
possible spanwise or separated flow effects on the outboard control surface. 
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as -^^scussed previously for the steatly-oressure data. The experimental and the 
unsteady-pressure magnitude data show good correlations In trends forward of 
each control -surface hinge line. Overall the deviation In this region for 
magnitudes Is more significant for the Inboard control surface. The Doublet 
Lattice phase angles were more negative (lagging) than the experimental data 
toward the leading edge. There was good correlation between theoretical and 
experimental phase angles In trend, with the best conyjarlson at the In-phase 
point ■ 0) of approximately 60 percent chord. There Is a noticeable 
deviation between the theoretical and experimental phase angles at the same 
ti^alllng edge locations of the outboard control surface which have shown 
possible spanwise or separated flow effects In the pressure magnitude data. 

CONCLUUiNG HLMARKS 

This paper presents comparisons of theoretical and experimental steady- am1 
unsteady-pressure distributions on a hlgh-aspect-ratio supercritical wing nwlel 
at M ■ 0.60. The theoretical calculations were performed using the unsteady 
aerodynamic 1 Ift^ng-surface method of Doublet Lattlc:. The trends of the 
Doublet Lattice data show an overall good comparison to both the steady and 
unsteady exp®» .inental aerodynamics. Comparisons for the steady data Include 
chordwise Incremental lifting-surface pressure distributions per Incremental 
angle of attacx and Incremental control -surface defle.'t i(.n, spanwise Incremental 
lift and moiitent distributions per Incremental angle of attack, spanwise 
distributions of local aerodynamic center locations, and spanwise Incremental 
lift distributions per Incremental Inboard and outboard control -surface 
deflection. The theoretical and experimental 1 Ifting-surface unsteady-pressure 
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coefficient magnitude and phase angle comparisons are shorn, for both Inboard and 
outboard oscillating control surfaces. 

Although, In general, the experimental and theoretical data show reasonably 
good agreement, there are some significant areas of deviation. These 
differences, which could possibly be adjusted with empirically-developed 
correction factors, are noted as follows: 

1. The Doublet Lattice data underpredicts experimental chordwise steady-* 
pressure distributions toward the leading edge and overpredi ts values toward 
the trailing edge. This deviation Is typical of thickness effects (viscosity) 
as discussed in ref. 8. 

2. The most significant deviation Is the Doublet Lattice overprediction of 
experimental chordwise pressure distributions on the control surfaces. The 
difference In both steady- and unsteady-pressure magnitudes Is more pronounced 
on the outboard control surface. 

3. For the unsteady aerodynamics, the Doublet Lattice phase angle values 
are consistently more negative (lagging) than the experimental values toward the 
leading edge. 

4. The experimental chordwise pressure distributions toward the trailing 
edge of the outboard control surface exhibited an effect possibly due to 
spanwise or separated flow. This Is Illustrated by a measured pressure rise at 
the trailing edge chordwise station x/c ■ 0.95 after a minimum pressure was 
reached forward of this location at x/c » 0.90. 
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TABLE I. 


- SUMMARY OF EXPERIMENTAL STEAOY-PRESSURE TEST CONDITIONS AT M 


Test Point No. 

a 

5 

(See Ref. 7) 

deg 

deg 


Inboard Control Surface 


265 

0 

+ 6 

266 



+ 4 

267 



+ 2 

268 



- 2 

269 



- 4 

270 



- 6 

271 



0 

186 

2.8^ 


+ 6 

187 



+ 4 

188 



+ 2 

190 



- 2 

191 



- 4 

192 



- 6 

193 



0 


Outboard Control Surface 


218 

219 

222 

223 

224 

225 

226 


+ 4 
+ 2 
- 2 

- 4 

- 6 
0 


174 
161 

175 
178 
162 
184 


2.85 


+ 

+ 

+ 


185 


6 

4 

2 

2 

4 

6 

0 
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T<VBLE II. 


- SUMMARY OF EXPERIMENTAL UNSTEADY-PRESSURE TEST CONDITIONS 


AT M • 0.60 


Test Point No. 

a 

5 

f 

k 

(See Ref. 7) 

deg 

deg 

Hz 

Inboard Control 

Surface 


121 

0 

± i 


5 

0.135 

124 

1 



10 

0.272 

125 

1 



15 

0.408 

145 

2.85 



5 

0.136 

144 

2.85 



10 

0.272 

143 

2.84 



15 

0.408 

j Outboard Control Surface 

130 

0 

t ^ 


5 

0.136 

129 

1 



10 

0.273 

128 

I 



15 

0.410 

139 

2.84 



5 

0.136 

140 

2.84 



10 

0.272 

142 

2.84 



15 

0.408 
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Orlfico Somlipon 

Stotiont 

lew No. 


1 

0.19 

2 

0.23 

3 

0.25 

4 

0.33 

5 

0.S1 

6 

0.71 

7 

0.75 

8 

0.51 

9 

0.92 



AR = 10.76 

S=1.94 m2 


% 


c,^=.4M 


Inboard Control Surfoco 


Outboard Control Surfoco 


2.2S6 


Sketch of wing planform geometry and orifice semispan stations 
Linear dimensions in meters. 






Figure 2.- Sketch of supercritical airfoil for 3-dimensional wind-tunnel 
model and plot of local chord and thickness variation along 
semi span. 
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Figure 3.- Sketch of Doublet Lattice aerodynamic model. 




EXPERIMENTAL 


OF POCU s 


IP^ 


O Test Points 185 & 226 

□ Test Points 193 & 271 



(a) At semi span location n * 0.19. 

Figure 4. - Chordwise incremental lifting-surface steady-pressure distribution 
for an incremental angle of attack. 
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EXPERIMENTAL 

O Test Points 185 & 226 

□ Test Points 193 & 271 


— + THEORETICAL 
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EXPERIMENTAL 


O Test Po1r*s 185 & 226 

□ Test Points 193 & 271 



0 .2 .k ,6 .8 1.0 

x/c 


(e) At semi span location n ■ 0.51. 
Figure 4,- Continued. 
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(f) At semispan location n ■ 0.71. 

Figure 4.- Continued. 




27 



OF fug:i « G'jvY 

EXPERIMENTAL 

O Test Points 185 A 226 

□ Test Points 193 A 271 


+ 


THEORETICAL 



OF POOH QUALITY 


EXPERIMENTAL 

O Test Points 18S « 226 

□ Test Points 193 I 271 



0 .2 .4 .6 .8 1.0 

x/c 

(h) At se«n1span location n * 0.81. 

Figure 4.- Continued. 
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(a) For A6 = ± 2° and o = 0°. 

Figure 7. - Chordwise incremental lifting-surface steacly-pressure distribution 
for an incremental inboard control -surface deflection at semispan 
station n = 0.19. 
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Figure 7.- Continued. 
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(f) For A6 ■ * 6® and O ■ 2.84®. 
Figure 7.- Concluded. 
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(a) For a 5 * t 2® and a * 0®. 

Figure 8. - Chordwise Incremental lifting-surface steady-pressure distribution 
for an incremental outboard control -surface deflection at semi span 
station n » 0.71. 
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(e) For A6 ■ t 4* and a • 2.85*. 
Figure 3.- Continued. 
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Figure 8.- Concluded. 
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(a) For 6 • t 4* and k • 0.14 (f • 5Hi). 

Figure 11. - Magnitude and phase angle of chordelse lifting-surface 

unsteady-pressure distribution at senispan station n ■ 0.71, due 
to oscillating outboard control surface. 
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